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Abstract: This paper presents the derivation of a new algorithm for multiplying of two Kaluza numbers. 

Performing this operation directly requires 1024 real multiplications and 992 real additions. The 
proposed algorithm can compute the same result with only 512 real multiplications and 576 real 
additions. The derivation of our algorithm is based on utilizing the fact that multiplication of 
two Kaluza numbers can be expressed as a matrixvector product. The matrix multiplicand that 
participates in the product calculating has unique structural properties. Namely exploitation of 
these specific properties leads to significant reducing of the complexity of Kaluza numbers multi¬ 
plication. 
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1. Introduction 

Today hypercomplex numbers [[Tl are used in various fields of data processing including 
digital signal and image processing, machine graphics, telecommunications and cryptography 
[S 111 111151 IB [71 [HIlKIOl. The most popular are quaternions, octonions and sedenions [[H. Per¬ 
haps the less popular are the Pauli, Dirac and Kaluza numbers. This numbers are mostly used 
in solving different physical problems in electrodynamics, field theory, etc. Somehow or other, 
hypercomplex number systems and their applications to data processing are beautiful enough to 
be worth studying simply for the pleasure of it. 
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The multiplication of two hypercomplex numbers is often one of the more time consuming 
operations. The reason for this is, because the addition of iV-dimensional hypercomplex numbers 
requires N real additions, while the multiplication of these numbers already requires N{N — 1) 
real additions and real multiplication. It is easy to see that the increasing of dimensions of 
hypercomplex number increases the computational complexity of the multiplication. Therefore, 
reducing the computational complexity of the multiplication of hypercomplex numbers is an im¬ 
portant scientific and engineering problem. 

Efficient algorithms for the multiplication of various hypercomplex numbers already exist 
[fT^ [T3l [T4l [T5l [m [nl [TSl [T9l [201 Ell E2l EH Ell. No such algorithms for the multiplication 
of Kaluza numbers have been proposed. The aim of the present paper is to suggest an efficient 
algorithm for this purpose. 

2. Preliminary Remarks 

A Kaluza number is defined as follows: 

31 

d do T ^ ^ dfiin, (2.1) 

r7.=l 

where {do} and {d„}, n = 1,..., 31 are real numbers, and {in}, n = 1,..., 31 are the imaginary 

units. The results of all possible products of Kaluza numbers imaginary units can be summarized 
in the following table: 


Table 1. North-west quadrant of the table for Kaluza numbers imaginary units multiplication 


X 

1 

61 

62 

63 

64 

65 

66 

67 

68 

eg 

610 

611 

612 

613 

614 

615 

1 

1 

61 

62 

63 

64 

65 

66 

67 

68 

eg 

610 

611 

612 

613 

614 

615 

ei 

6l 

1 

66 

67 

68 

69 

62 

63 

64 

65 

616 

617 

618 

619 

620 

621 

62 

62 

-66 

1 

610 

611 

612 

-6l 

-616 

-617 

-618 

63 

64 

65 

622 

623 

624 

63 

63 

-67 

-6l0 

-1 

613 

614 

6i6 

6l 

-619 

— 620 

62 

— 622 

— 623 

-64 

-65 

625 

64 

64 

-68 

-611 

-613 

-1 

615 

617 

619 

61 

— 621 

622 

62 

— 624 

63 

— 625 

-65 

65 

65 

-69 

— 6i2 

— 6i4 

— 6i5 

-1 

618 

620 

621 

61 

623 

624 

62 

625 

63 

64 

66 

66 

-62 

61 

616 

617 

618 

-1 

-610 

-611 

— 612 

67 

68 

eg 

626 

627 

628 

67 

67 

-63 

-616 

-6l 

619 

620 

6l0 

1 

-613 

— 614 

66 

— 626 

— 627 

-68 

-69 

629 

68 

68 

-64 

-617 

-619 

-6l 

621 

6ll 

613 

1 

— 615 

626 

66 

— 628 

67 

— 629 

-69 

69 

69 

-65 

-618 

— 620 

— 621 

-6l 

612 

614 

615 

1 

627 

628 

66 

629 

67 

68 

6l0 

6l0 

616 

-63 

-62 

622 

623 

-67 

-66 

626 

627 

1 

-613 

— 614 

-611 

— 612 

630 

6ll 

6ll 

617 

-64 

— 622 

-62 

624 

-68 

— 626 

-66 

628 

613 

1 

— 615 

610 

-630 

— 612 

6i2 

6i2 

618 

-65 

— 623 

— 624 

-62 

-69 

— 627 

— 628 

-66 

614 

615 

1 

630 

610 

611 

6i3 

6i3 

619 

622 

64 

-63 

625 

626 

68 

-67 

629 

611 

-610 

630 

-1 

615 

— 614 

6i4 

6i4 

620 

623 

65 

— 625 

-63 

627 

69 

— 629 

-67 

612 

-630 

-610 

— 615 

-1 

613 

615 

615 

621 

624 

625 

65 

-64 

628 

629 

eg 

-68 

630 

612 

-611 

614 

-613 

-1 
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Table 2. North-east quadrant of the table for Kaluza numbers imaginary units multiplieation 


X 

6i6 

617 

6i8 

619 

620 

621 

622 

623 

624 

625 

626 

627 

628 

629 

630 

631 

1 

6i6 

617 

6i8 

619 

620 

621 

622 

623 

624 

625 

626 

627 

628 

629 

630 

631 

ei 

6l0 

611 

612 

613 

614 

615 

626 

627 

628 

629 

622 

623 

624 

625 

631 

630 

62 

-67 

-68 

-69 

— 626 

—627 

— 628 

— 613 

614 

615 

630 

— 619 

— 620 

—621 

— 631 

625 

—629 

63 

-66 

626 

627 

68 

69 

— 629 

611 

612 

-630 

-615 

-617 

-618 

631 

621 

624 

—628 

64 

— 626 

-66 

628 

-67 

629 

69 

— 610 

630 

612 

614 

6i6 

— 631 

—618 

— 620 

— 623 

627 

65 

— 627 

—628 

-66 

— 629 

-67 

-68 

-630 

—610 

-611 

—613 

631 

616 

617 

619 

622 

—626 

66 

-63 

-64 

-65 

— 622 

—623 

— 624 

619 

620 

621 

631 

— 613 

— 614 

-615 

— 630 

629 

—625 

67 

-62 

622 

623 

64 

65 

-625 

617 

6i8 

— 631 

-621 

-611 

-612 

630 

615 

628 

—624 

68 

— 622 

-62 

624 

-63 

625 

65 

— 6i6 

631 

6i8 

620 

610 

-630 

-612 

— 614 

— 627 

623 

69 

— 623 

—624 

-62 

— 625 

-63 

-64 

— 631 

—616 

-617 

—619 

630 

610 

611 

613 

626 

—622 

610 

61 

—619 

—620 

-617 

—618 

631 

64 

65 

-625 

—624 

68 

69 

—629 

— 628 

615 

621 

611 

6i9 

61 

-621 

616 

—631 

— 6i8 

-63 

625 

65 

623 

-67 

629 

69 

627 

— 614 

—620 

612 

620 

621 

61 

631 

616 

617 

-625 

-63 

-64 

—622 

—629 

-67 

-68 

— 626 

613 

619 

613 

617 

—616 

631 

-61 

621 

— 620 

-62 

624 

— 623 

-65 

-66 

628 

—627 

-69 

— 612 

—618 

614 

6i8 

—631 

—616 

-621 

-61 

619 

— 624 

-62 

622 

64 

— 628 

-66 

626 

68 

611 

617 

615 

631 

618 

-617 

620 

—619 

-61 

623 

—622 

-62 

-63 

627 

— 626 

-66 

-67 

— 610 

—616 


Table 3. South-west quadrant of the table for Kaluza numbers imaginary units multiplieation 


X 

1 

61 

62 

63 

64 

65 

66 

67 

68 

69 

610 

611 

612 

613 

614 

ei5 

616 

616 

610 

-67 

-66 

626 

627 

-63 

-62 

622 

623 

61 

— 619 

— 620 

-617 

— 618 

631 

617 

617 

611 

-68 

— 626 

-66 

628 

-64 

— 622 

-62 

624 

619 

61 

— 621 

616 

— 631 

— 618 

618 

618 

612 

-69 

— 627 

— 628 

-66 

— 65 

— 623 

— 624 

-62 

620 

621 

61 

631 

616 

617 

619 

619 

613 

626 

68 

-67 

629 

622 

64 

-63 

625 

617 

— 616 

631 

-61 

621 

— 620 

620 

620 

614 

627 

69 

— 629 

-67 

623 

65 

-625 

-63 

618 

— 631 

— 616 

— 621 

-61 

619 

621 

621 

615 

628 

629 

69 

-68 

624 

625 

65 

-64 

631 

618 

-617 

620 

— 619 

-61 

622 

622 

— 626 

613 

611 

— 610 

630 

— 619 

-617 

616 

— 631 

64 

-63 

625 

-62 

624 

— 623 

623 

623 

— 627 

614 

612 

-630 

— 610 

— 620 

— 618 

631 

616 

65 

-625 

-63 

— 624 

-62 

622 

624 

624 

— 628 

615 

630 

612 

-611 

— 621 

— 631 

— 618 

617 

625 

65 

-64 

623 

— 622 

-62 

625 

625 

— 629 

-630 

— 615 

614 

— 613 

631 

621 

— 620 

619 

624 

—623 

622 

-65 

64 

-63 

626 

626 

— 622 

619 

617 

— 616 

631 

— 613 

-611 

610 

-630 

68 

-67 

629 

-66 

628 

— 627 

627 

627 

— 623 

620 

618 

— 631 

— 616 

— 614 

— 612 

630 

610 

69 

—629 

-67 

— 628 

-66 

626 

628 

628 

— 624 

621 

631 

618 

-617 

— 615 

-630 

— 612 

611 

629 

eg 

-68 

627 

— 626 

-66 

629 

629 

-625 

— 631 

— 621 

620 

— 619 

630 

615 

— 614 

613 

628 

— 627 

626 

-69 

68 

-67 

630 

630 

631 

-625 

— 624 

623 

— 622 

— 629 

— 628 

627 

— 626 

615 

— 614 

613 

— 612 

611 

— 610 

631 

631 

630 

— 629 

— 628 

627 

— 626 

-625 

— 624 

623 

— 622 

621 

— 620 

619 

— 618 

617 

— 616 
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Table 4. South-east quadrant of the table for Kaluza numbers imaginary units multiplication 


X 

616 

617 

618 

619 

620 

621 

622 

623 

624 

625 

626 

627 

628 

629 

630 

631 

ei6 

1 

— 613 

—614 

-611 

—612 

630 

68 

69 

— 629 

— 628 

64 

65 

—625 

— 624 

621 

615 

ei 7 

613 

1 

-615 

610 

-630 

-612 

-67 

629 

69 

627 

-63 

625 

65 

623 

— 620 

—614 

ei8 

614 

615 

1 

630 

610 

611 

— 629 

-67 

-68 

— 626 

—625 

-63 

-64 

— 622 

619 

613 

ei 9 

611 

— 610 

630 

-1 

615 

— 614 

-66 

628 

— 627 

-69 

-62 

624 

—623 

-65 

— 618 

—612 

620 

612 

-630 

—610 

-615 

-1 

613 

— 628 

-66 

626 

68 

—624 

-62 

622 

64 

617 

611 

621 

630 

612 

-611 

614 

—613 

-1 

627 

—626 

-66 

-67 

623 

— 622 

-62 

-63 

— 616 

—610 

622 

-68 

67 

—629 

66 

—628 

627 

-1 

615 

— 614 

-612 

61 

-621 

620 

618 

-65 

69 

623 

-69 

629 

67 

628 

66 

— 626 

-615 

-1 

613 

611 

621 

61 

—619 

-617 

64 

-68 

624 

—629 

-69 

68 

— 627 

626 

66 

614 

—613 

-1 

— 610 

—620 

619 

61 

616 

-63 

67 

625 

— 628 

627 

—626 

69 

-68 

67 

612 

-611 

610 

1 

—618 

617 

—616 

-61 

-62 

66 

626 

-64 

63 

—625 

62 

—624 

623 

-61 

621 

— 620 

— 618 

1 

-615 

614 

612 

-69 

65 

627 

-65 

625 

63 

624 

62 

— 622 

— 621 

-61 

619 

617 

615 

1 

—613 

-611 

68 

-64 

628 

-625 

-65 

64 

— 623 

622 

62 

620 

—619 

-61 

— 616 

—614 

613 

1 

610 

-67 

63 

629 

— 624 

623 

—622 

65 

-64 

63 

618 

-617 

616 

61 

-612 

611 

—610 

-1 

-66 

62 

630 

621 

— 620 

619 

— 618 

617 

— 616 

65 

-64 

63 

62 

69 

-68 

67 

66 

-1 

-61 

631 

615 

— 614 

613 

-612 

611 

— 610 

69 

-68 

67 

66 

65 

-64 

63 

62 

-61 

-1 
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Suppose we want to compute the product of two Kaluza numbers. 


31 

d = = do + dnim 

n=l 


where 


31 31 

d^^^ = no + ^ anin, d^‘^^ = ^ Kin- 

n=l n=l 


( 2 . 2 ) 


The operation of multiplication of Kaluza numbers can be represented more compactly in the 
form of vector-matrix product: 


Y32XI — B32X32XI, 


(2.3) 


where 


X32XI — [flO) Ol, • • • , asi]"*", 
Y32XI = [do, di,..., dis]"*", 


B32 — 


■•^le ■'^16 

(0,1) g(l,l) 


B 


16 


16 


bo 

bi 

bo 

-bo 

-bi 

—bo 

-bo 

67 

bs 

bo 

bio 

bn 

bio 

—bio 

—bii 

—bio 

bi 

bo 

-bo 

67 

bs 

bo 

bo 

-bo 

-bi 

—bo 

bio 

bn 

bis 

—bio 

—boo 

—boi 

&2 

bo 

bo 

bio 

bn 

bio 

-bi 

—bio 

—bn 

—bis 

-bo 

-bi 

—bo 

—boo 

—boo 

—boi 

bo 

67 

bio 

bo 

bio 

bii 

—bio 

-bi 

—bio 

—boo 

-bo 

—boo 

—boo 

-bi 

—bo 

-boo 

64 

bs 

bn 

—bio 

bo 

bio 

—bn 

bio 

-bi 

—boi 

boo 

-bo 

—boi 

bo 

boo 

—bo 

bo 

bo 

bio 

—bii 

—bio 

bo 

—bis 

boo 

boi 

-bi 

boo 

boi 

-bo 

-boo 

bo 

bi 

bo 

bo 

-bi 

—bie 

—bn 

—bis 

bo 

bio 

bn 

bio 

—67 

-bs 

-bo 

—boo 

—bo 7 

—bos 

67 

bo 

—biQ 

-bi 

—bio 

—boo 

bio 

bo 

bio 

bii 

—bo 

—boo 

—bo 7 

-bs 

-bo 

—boo 

bs 

bi 

—bn 

bio 

-bi 

—boi 

bn 

—bio 

bo 

bio 

boo 

—bo 

—bos 

67 

boo 

-bo 

bo 

bo 

—bis 

boo 

boi 

-bi 

bio 

—bii 

—bi 5 

bo 

bo 7 

bos 

-bo 

—boo 

67 

bs 

bw 

bio 

bo 

-bo 

—boo 

—boo 

—67 

bo 

boo 

bo 7 

bo 

bio 

bii 

-bn 

—bio 

—boo 

bn 

bn 

bi 

boo 

-bo 

—boi 

-bs 

—boo 

bo 

bos 

—bio 

bo 

bio 

bio 

boo 

—bio 

b \2 

bio 

bo 

boo 

boi 

-bo 

-bo 

— bo 7 

—bos 

bo 

—bii 

—bio 

bo 

—boo 

bio 

bn 

bi 3 

bio 

boo 

bi 

-bo 

-boo 

—boo 

-bs 

67 

boo 

-bn 

bio 

boo 

bo 

bio 

—bii 

bi 4 

boo 

boo 

bo 

boo 

-bo 

—bo 7 

-bo 

—boo 

67 

—bio 

—boo 

bio 

—bio 

bo 

bio 

bi 5 

boi 

boi 

-boo 

bo 

-bi 

—bos 

boo 

-bo 

bs 

boo 

—bio 

bn 

bii 

—bio 

bo 
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bo 

bi3 

614 

-bn 

—bio 

—boo 

bs 

bo 

boo 

00 

1 

-b4 

—bs 

-bos 

boA 

—boi 

-bis 

—bi3 

bo 

bis 

bio 

boo 

—bio 

—67 

—boo 

bo 

bo7 

bo 

bos 

—bs 

—boo 

boo 

biA 

— &14 

-bis 

bo 

—boo 

bio 

bn 

boo 

— &7 

—bs 

—boo 

-bos 

bo 

bA 

boo 

—bio 

—bio 

-bn 

bio 

boo 

bo 

bis 

—bi4 

-bo 

—bos 

bo7 

bo 

bo 

boA 

—boo 

-bs 

bis 

bio 

—b\2 

—boo 

bio 

-bis 

bo 

bio 

bos 

—bo 

-boo 

—bs 

—bo4 

bo 

boo 

bA 

—bi7 

-bn 

boo 

—bio 

bn 

bi4 

—bio 

bo 

—bo7 

boo 

-bo 

67 

boo 

—boo 

bo 

—bo 

biG 

bio 

bs 

—67 

—boo 

bo 

bos 

— bo7 

bo 

bis 

— &14 

bio 

-bi 

—boi 

boo 

—bis 

-bs 

-bo 

bo 

boo 

—67 

—bos 

bo 

boo 

-bis 

bo 

bio 

-bn 

boi 

-bi 

—bio 

bi7 

bA 

bs 

—boo 

bo 

-bs 

bo7 

—boo 

bo 

bi4 

—bio 

bo 

bio 

—boo 

bio 

-bi 

—bio 

-bo 

—67 

—bos 

bo7 

—boo 

bo 

-bs 

67 

bio 

-bn 

bio 

bo 

—bis 

bi7 

—bio 

-bi 

-bo 

-bo 

bi 

—bo 

-bos 

bo 

bo4 

—boo 

-bi 

—boi 

boo 

“^18 

bo 

bis 

—biA 

bio 

—bo 

-bs 

bs 

bos 

-bo 

—bo4 

bo 

boo 

boi 

-bi 

—bio 

&17 

-bis 

bo 

bio 

-bn 

bs 

b4 

—bos 

bs 

-bi 

boo 

—boo 

bo 

—boo 

biQ 

-bi 

— ^16 

bi4 

—bio 

bo 

bio 

—67 

-bo 

—boi 

boo 

—boo 

bs 

-64 

bo 

—bis 

bi7 

—bio 

-bi 

bio 

-bn 

bio 

bo 

-bo 

-bo 

boi 

—boo 

bio 

bis 

—bi7 

bie 

bs 

-64 

bo 

—bo 

—bo 

bs 

—67 

bo 

bo 

bi 

bis 

— &14 

bio 

bio 

-bn 

bio 

—bo 

bs 

— &7 

bo 

bs 

-bA 

bo 

-bo 

bi 

bo 


3. Synthesis of a rationalized algorithm for computing Kaluza numbers 
product 

The direct multiplication of two Kaluza requires 1024 real multiplications and 992 real addi¬ 
tions. We shall present the algorithm, which reduce arithmetical complexity to 512 real multipli¬ 
cations and 576 real additions. 

At first, we rearrange the rows of the matrix in the following order {1, 2, 3, 7, 5, 9, 4, 8, 6, 
10, 11, 17, 13, 19, 15, 21, 12, 18, 14, 20, 16, 22, 23, 27, 25, 29, 24, 28, 26, 30, 31, 32}. Next, 
we rearrange the columns of obtained matrix in the same manner. As a result, we obtain the 
following matrix: 

Then we can rewrite expression (2.3) in following form: 


Y32XI — P32B32P32X32X1, 


(3.1) 


32 — 


p{0i0) p(liO) 
■'^16 

Td{0i1) T3(1i1) 


16 


16 


where 



1 0 0 0 0 0 

0 1 0 0 0 0 

0 0 1 0 0 0 

0 0 0 0 0 0 

0 0 0 0 1 0 

0 0 0 0 0 0 

0 0 0 1 0 0 

p{o,o) _ 0 0 0 0 0 0 

- 0 0 0 0 0 1 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

_ 0 0 0 0 0 0 

'000000 
0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

p{o,i) _ 0 0 0 0 0 0 

16“ 000000 
0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

'000000 
0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

p{i,o) _ 0 0 0 0 0 0 

16“ 000000 
0 0 0 0 0 0 

0 0 0 0 0 0 

1 0 0 0 0 0 

0 0 0 0 0 0 

0 0 1 0 0 0 

0 0 0 0 0 0 

0 0 0 0 1 0 


0000000000 
0000000000 
0000000000 
1000000000 
0000000000 
0010000000 
0000000000 
0100000000 
0000000000 ’ 
0001000000 
0000100000 
0000000000 
0000001000 
0000000000 
0000000010 
0000000000 _ 

0000010000 ' 

0000000000 

0000000100 

0000000000 
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If we interpret the resulting matrix as a bloek-matrix, it is easy to see that eaeh bloek of this 
matrix (or more preeisely, eaeh (2 x 2 )-submatrix) is bisymmetrie, i.e. has the properties of pair 
wise symmetry about both of its main diagonals. There is an effeetive method of faetorization 
of this type matriees, whieh allows during the ealeulation of the veetor-matrix produets halve the 
number of multiplieations, that is to perform only two, not four multiplieations at the expense of 
triple inereasing the number of additions (from two to six) [l25ll2^ : 


a 

b ■ 


Xq 


■ 1 

1 


|(ct -|- 6) 

0 


■ 1 

1 


Xq 

b 

a 


. . 


1 

-1 


0 

|(a-6) _ 


1 

-1 


. . 


where 


= Ha is the (2 X 2) Hadamard matrix. 


1 1 
1 -1 

Fig. 1. shows a data flow diagram of the rationalized algorithm for (2 x 2) bisymmetrie 
matrix-veetor multiplieation in aeeording to (3.2). The eireles in this figure show the opera¬ 
tion of multiplieation by a value inseribed inside a eirele. In turn, the reetangles indieate the 
matrix-veetor multiplieations with matriees inseribed inside reetangles. In this paper the data 
flow diagrams are oriented from left to right. Straight lines in the figure denote the operation of 
data transfer. We use the solid lines without any arrows, so as not to elutter up the presented 
diagrams. 

Now it is easy to show that using the above method of (2 x 2)-bisymmetrie matrix faetorization 
we ean eonstruet an effieient algorithm for multiplying Kaluza numbers. In sueh algorithm the 
number of real multiplieations will be redueed by half eompared with naive method of ealeu- 
lations. However, if we direetly use an expansion (3.2) for eaeh of the submatriees of matrix 
B 32 , the number of real additions will inereases dramatieally. In this ease, the total value of 
eomputation eomplexity beeomes even greater than in the ease of naive method of eomputations. 
However, if we utilize the faet that the multiplieation of the (2 x 2) -Hadamard matrix on the 
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corresponding subvector (the fig. 1 is positioned before of multiplications) is a common opera¬ 
tion for all submatrices disposed in a same vertical of matrix B 32 , and a similar operation (in the 
fig.l, it follows by multiplications) is a common operation for all submatrices disposed in a same 
horizontal of matrix B 32 , the number of real additions can be significantly reduced. Consider the 
synthesis of an efficient algorithm for multiplying Kaluza numbers in more detail. 


Xq 

Xi 


H, 


Ki>& 


.yo 


So = l{a + b) Si = |(a - 6) 


Figure 1. Data flow diagram of the rationalized algorithm for (2 x 2)-bisymmetric matrix-vector 

multiplication in according to (3.2) 


Let us first introduce some matrices 

W 16 = (I 16 ® H 2 ), P512x32 = I 16 ® (Il6xl ® I 2 ), 

15 

D 512 = 032 ^ A 32 X 512 = (Il6xl C) I 32 ), 

k=0 

where Ia? is an identity N x N matrix, signs and “©” denote tensor product and direct sum 
of two matrices respectively, and lArxM is an integer matrix consisting of all Is [26], 

= diag{s^Q\ 4?), A; = 0,1,..., 15, 

where 
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*20 ~ 

C3O) 

*21 ~ 

C3I, 

*22 ~ 

C4, 

*23 ~ 

CS) 

- ■ 
*24 ~ 

“C28, 

- ■ 

*25 ~ 

“C29, 

*26 ~ 

CSi 

*27 ~ 

C9, 

- ■ 
*28 ~ 

“C24, 

- ■ 

*29 ~ 

“C25, 

*30 ~ 

C20; 

*31 “ 

C21, 


*0 — 

C12, 

*1 — 

Cl3) 

*2 — 

C8, 

*3 — 

~Cg, 

- ■ 
*4 — 

“C24, 

- ■ 

*5 — 

“C25, 

- ■ 
*6 — 

“C26, 

- ■ 
*7 — 

“C27 

*8 ~ 

“C2, 

*9 ~ 

“C3, 

c(6) - 
*10 ~ 

Cl4, 

*11 ~ 

Cl5, 

*12 ~ 

■ Co, 

*13 ~ 

^ Cl, 

*14 ~ 

ClO, 

*15 “ 

Cll, 

*16 ~ 

C2O) 

*17 ~ 

C2I) 

*18 ~ 

C30, 

*19 ~ 

C3I, 

*20 ~ 

C16) 

*21 ~ 

Cl7) 

- ■ 

*22 ~ 

“C28, 

- ■ 
*23 ~ 

“C29 

5(6) _ 
*24 ~ 

— C4, 

*25 ~ 

“C5, 

.(6) - 
*26 ~ 

“Cg, 

*27 ~ 

-C7, 

- - 
*28 ~ 

“C22, 

- - 

*29 ~ 

“C23, 

*30 ~ 

C18, 

*31 ~ 

Cl9, 
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- - 
*0 — 

“Ci 4 , 

- - 
*1 — 

“Ci 5 , 

- - 
*2 — 

“C26, 

- - 
*3 — 

“C27 

64 — 

C28) 

*5 — 

C29, 

*6 — 

^8; 

*7 — 

— Cg, 

II 

Cg, 

*9 — 

C7, 

- - 
*10 ~ 

“C12, 

- - 
*11 — 

-Cl 3 

*12 — 

ClO) 

s^^'> - 
*13 ~ 

Cll, 

*14 — 

Co, 

*15 — 

Cl, 

*16 ~ 

C30) 

J 7 ) _ 

*17 — 

C3I, 

*18 ~ 

C20, 

*19 ~ 

C2I, 

- - 
*20 ~ 

“C18, 

- - 
*21 — 

“Ci 9 , 

*22 — 

C24, 

*23 ~ 

C25, 

- - 
*24 — 

-C22, 

- - 
*25 ~ 

“C23, 

*26 ~ 

C2, 

*27 ~ 

C3, 

*28 ~ 

C4, 

*29 — 

^ 5 , 

- - 
*30 — 

~Cig, 

- - 
*31 ~ 

“Ci 7 


(8) 

So = CI6, 

*f) = 

Cl 7 , 

*2 ~ 

— C4, 

4 *’ = 

C5, 

J8) _ „ 

S4 — ~C2, 

*f = 

-C3, 

- ■ 

*6 ~ 

-C22, 

4*' = - 

-C23 

(8) 

Ss = C24, 

*9 ~ 

C25, 

*10 ~ 

C18, 

s‘®> - 
*11 ~ 

Cl 9 , 

(8) 

*12 = -C20, 

- - 
*13 ~ 

“C21, 

- ■ 

*14 ~ 

“C30, 

- ■ 

*15 ~ 

-C31 

s 

II 

s® - 

*17 ~ 

Cl, 

*18 ~ 

ClO, 

*19 ~ 

Cll, 

(8) 

*y = -C12, 

- - 
*21 ~ 

“Ci 3 , 

*22 ~ 

Cg, 

*23 ~ 

-C7, 

J8) _ „ 

*24 — Cs, 

*25 ~ 

C9, 

*26 ~ 

C28, 

*27 ~ 

C29, 

(8) 

*28 — ** 26 , 

*29 ~ 

C27, 

- ■ 

*30 ~ 

“Ci 4 , 

- ■ 

*31 ~ 

“C15 


( 9 ) 

*0 = -C18, 

( 9 ) 

*1 = -Cl 9 , 

*2 — 

— C22, 

*3 — 

— C23 

( 9 ) 

*4 = Cg, 

( 9 ) 

*5 = C 7 , 

*6 — 

— C4, 

*7 — 

C5, 

( 9 ) 

*8 = -C 28 , 

( 9 ) 

*9 = -C29, 

*10 ~ 

“Cig, 

*11 ~ 

“Ci 7 

( 9 ) 

*12 — “C30, 

( 9 ) 

*13 = -C3I, 

*14 ~ 

“C20, 

*15 ~ 

-C21 

( 9 ) 

*lg = ClO, 

J 9 ) _ „ 

*17 ~ Cll, 

*18 “ 

= Co, 

*19 “ 

= Cl, 

( 9 ) 

*20 — ** 14 , 

( 9 ) 

*21 — Cl 5 , 

*22 “ 

= C2, 

*23 “ 

= C3, 

( 9 ) 

*24 — C2g, 

( 9 ) 

*25 ~ * 727 , 

*26 ~ 

~C 24 , 

*27 ~ 

~C 25 

J 9 ) _ „ 

*28 — Cs, 

J 9 ) _ „ 

*29 — 4 - 9 , 

*30 “ 

' C12, 

s® - 

*31 “ 

' Ci 3 , 
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JIO) _ 
*0 ~ ■ 

“C20, 

Jio) _ 
*1 — ■ 

“C21, 

JIO) _ 
*2 ~ 

~C24, 

Jio) . 

*3 

JIO) _ 

64 — 

-Cs, 

JIO) _ 

*5 “ 

C9, 

Jio) _ 
*6 ~ 

-C28, 

Jio) . 

*7 

JIO) _ 
*8 ~ 

C4, 

Jio) _ 

*9 ~ 

C5, 

Jio) _ 
*10 ~ 

“C30, 

Jio) . 
*11 ■ 

II 

C16, 

Jio) _ 

*13 ~ 

Cl7, 

JIO) _ 

*14 ~ 

C18, 

Jio) 

*15 

JIO) _ 
^16 ~ ■ 

“C12, 

Jio) _ 

*17 “ ■ 

“C13, 

Jio) _ 
*18 ~ 

~Ci4, 

Jio) . 

*19 ■ 

0 :. 

to 

II 

Co, 

Jio) _ 
*21 ~ 

Cl, 

Jio) _ 
*22 ~ 

“C26, 

Jio) . 

*23 ■ 

Jio) _ 

'^24 “ 

C2, 

Jio) _ 

*25 ~ 

C3, 

Jio) _ 
*26 ~ 

C22, 

„(10) 

*27 

0 :. 

to 

II 

Ce, 

Jio) _ 

*29 ~ 

C7, 

Jio) _ 

*30 ~ 

ClO, 

Jio) 

*31 


*0 ~ 

-C23, 

=- 

-C22, 

5^“^ - 
*2 ~ 

“C19, 

- 

*3 

*4 — 

Cll, 

= 

ClO, 

*6 ~ 

“C17, 

- 

*7 

,(11) - . 
*8 — 

-C31, 

411) -. 

*9 — 

-C30, 

,(11) _ 
*10 ~ 

-C5, 

,(11) . 
*11 

o(ll) - . 
*12 — 

-C29, 

,(11) - . 
*13 “ 

-C28, 

,(11) _ 
*14 — 

— C25, 

,(11) _ 
*15 - 

*16 “ 

C7, 

*17 “ 

C6, 

*18 “ 

= C3, 

*19 

*20 ~ 

C27, 

*21 ~ 

C26, 

*22 “ 

= Cl, 

*23 

*24 ~ 

Cl5, 

*25 ~ 

Cl4, 

*26 ~ 

— C21, 

*27 " 

*28 ~ 

Cl3, 

*29 ~ 

C12, 

*30 “ 

= Cg, 

*31 

J12) _ 
*0 ~ ■ 

-C25, 

,(12) _ 
*1 — ■ 

-C24, 

,(12) _ 
*2 ~ 

— C21, 

,(12) _ 
*3 

J12) _ 

*4 — ■ 

-Ci3, 

,(12) _ 
*5 “ ■ 

-C12, 

,(12) _ 
*6 ~ 

— C31, 

,(12) _ 
*7 

J12) _ 
*8 ~ 

Cl7, 

,(12) _ 
*9 ~ 

C16, 

,(12) _ 
*10 ~ 

“C29, 

,(12) _ 
*11 " 

J12) _ 
*12 “ 

C5, 

,(12) _ 
*13 “ 

C4, 

,(12) _ 
*14 “ 

= C23, 

,(12) 

*15 

J12) _ 
*16 ~ 

“C9, 

,(12) _ 
*17 “ 

-C8, 

,(12) _ 
*18 ~ 

“C27, 

,(12) _ 
*19 " 

J12) _ 
*20 “ 

C3, 

,(12) _ 
*21 “ 

C2, 

,(12) _ 
*22 ~ 

“C15, 

,(12) _ 
*23 " 

J12) _ 

*24 “ 

Cl, 

,(12) _ 
*25 “ 

Co, 

,(12) _ 
*26 “ 

= C19, 

,(12) 

*27 

J12) _ 
*28 ~ 

ClO, 

,(12) _ 
*29 ~ 

Cll, 

,(12) _ 
*30 “ 

= C7, 

,(12) 

*31 


“C25, 
“C29, 
“C31, 
= Ci9, 
“Ci 5 , 

“C27, 

= C 23 , 

= Cii, 


“C 18 , 
“C16, 
-C4, 
— C24, 
= C2, 

= Co, 
—C20, 
= Cg, 


—C20, 
—C30, 
—C28, 
= C22, 
“C26, 

“Ci4, 

= C18, 

= Ce, 
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.(13) 

’0 ~ ■ 

“C27, 

„(13) _ 
*1 — ■ 

“C26, 

J13) _ 
*2 ~ ■ 

“Cl5, 

<,(13) _ 

*3 ~ ■ 

“Ci4, 

J13) _ 

*4 — 

C31j 

J13) _ 
*5 ~ 

C3O: 

J13) _ 
*6 ~ ■ 

“Cl3, 

<,(13) _ 
*7 ~ ■ 

“C12, 

J13) _ 
*8 ~ 

Cllj 

J13) _ 

*9 ~ 

ClOi 

J13) _ 
*10 ~ 

“C9j 

,(13) _ 
*11 ~ 

-C8, 

II 

’ C7, 

J13) _ 
*13 “ 


J13) _ 

*14 “ 

C3, 

,(13) _ 

*15 “ 

C2, 

J13) _ 
*16 ~ 

C29j 

J13) _ 

*17 ~ 

C 28 , 

J13) _ 
*18 ~ 

C25, 

,(13) _ 
*19 ~ 

C24, 

.(13) 

>20 ~ ■ 

-C23, 

J13) _ 
*21 “ ■ 

-C 22 , 

<,(13) _ 
*22 ~ 

C2I, 

,(13) _ 

*23 ~ 

C20, 

.(13) 

’24 ~ ■ 

“Cl9, 

„(13) _ 

*25 “ ■ 

“CI8, 

<,(13) _ 
*26 “ 

Cl, 

,(13) _ 
*27 “ 

Co, 

.(13) 

’28 ~ ■ 

-c-n, 

J13) _ 

*29 “ ■ 

“CI6, 

<,(13) _ 

*30 ~ 

-C5, 

,(13) _ 
*31 ~ 

-C4, 


,(14) _ 
*0 ~ 

C29, 

,(14) _ 
*1 ~ 

C28, 

,(14) _ 
*2 ~ ■ 

-C3I, 

,(14) _ 
*3 ~ 

-C30, 

II 

Cl5, 

,(14) _ 

*5 ~ 

Cl4, 

,(14) _ 
*6 ~ 

— C21, 

,(14) _ 

*7 ~ ■ 

-C20, 

,(14) 

8 ~ ■ 

-Cl9, 

,(14) _ 

*9 “ ■ 

-C18, 

,(14) _ 
*10 ~ 

C25, 

,(14) _ 
*11 ~ 

C24, 

II 

C23, 

,(14) _ 

*13 ~ 

C22, 

,(14) _ 

*14 “ 

= C5, 

,(14) _ 
*15 “ 

= C4, 

,(14) 

16 ~ ■ 

-C27, 

,(14) _ 

*17 “ ■ 

-C26, 

,(14) _ 
*18 ~ 

— C9, 

,(14) _ 
*19 ~ 

— C8, 

.(14) 

’20 ~ 

-C7, 

,(14) _ 
*21 “ 

Cg, 

,(14) _ 
*22 ~ 

Cl3, 

,(14) _ 
*23 ~ 

C12, 

,(14) _ 

*24 ~ 

Cll, 

,(14) _ 
*25 ~ 

Cio, 

,(14) _ 
*26 ~ 

-Ci7, 

,(14) _ 
*27 ~ ■ 

-C16, 

,(14) _ 
*28 “ 

’ Cl, 

,(14) _ 
*29 “ 

’ Co, 

,(14) _ 
*30 ~ 

-C3, 

,(14) _ 
*31 ~ 

-C2, 


(15) 

So = -C30, 

(15) 

*1 = -C3I, 

,(15) _ 
*2 ~ 

C28, 

(15) 

*3 = C29, 

(15) 

S4 = -C26, 

(15) 

*5 = -C27, 

,(15) _ 
*6 ~ 

C24, 

,(15) _ „ 

*7 — C25, 

(15) 

Ss = C22, 

(15) 

*9 = C23, 

,(15) _ 
*10 ~ ■ 

“C20, 

,(15) _ „ 

*11 — —C21, 

(15) 

Si2 = -C18, 

(15) 

*13 = -Cl9, 

,(15) _ 

*14 ~ ■ 

“C16, 

(15) 

*15 = -Cl7, 

(15) 

S16 = Cl4, 

,(15) _ „ 

*17 — C15, 

,(15) _ 
*18 ~ 

C12, 

(15) 

*19 = Cl3, 

(15) 

*20 — ClO, 

(15) 

*21 = Cll, 

,(15) _ 
*22 ~ 

-C8, 

,(15) _ „ 

*23 — C9, 

*24^^ = -C6, 

(15) 

*25 = -C7, 

,(15) _ 
*26 ~ 

C4, 

(15) 

*27 = C5, 

,(15) _ „ 

'^28 ~ ^2 5 

,(15) _ „ 

*29 ~ C3, 

,(15) _ 

*30 ~ 

Co, 

(15) 

*31 = Cl 


and 
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Co = 

K^o + bi ), 

Cl = 

|(6o 

- bi ), 

C2 = 

2(62 + bg ), 

C3 = 

1(62- 

- bg ). 

C4 = 

1(64 + bs ), 

C5 = 

1(64 ■ 

- bs ), 

Co = 

1(63 + 67), 

C7 = 

\{ b:i ■ 

- bj ). 

C8 = 

|(^5 + bg ), 

cg = 

1(^5 

- bg ), 

Cio = 

|(6io + 6ie), 

Cll = 

2(^10 

— big 

C12 = 

|(6i2 + bis ), 

Cl3 = 

|(6i2 

- bis ), 

Cl4 = 

2(614 + 620), 

Cl5 = 

2(^14 

— 620 

C16 = 

l(611 + 617), 

C17 = 

|(6ii 

— bn ), 

C18 = 

2(613 + 619), 

Cl9 = 

2(^13 

~ 619 

C20 = 

1(^15 + ^21), 

C21 = 

|(6 i 5 

— bgi ). 

C22 = 

1(622 + 620), 

C23 = 

2(^22 

— 626 

C24 = 

1(624 + bgs ), 

C25 = 

2(624 

— bgs ). 

C26 = 

1(623 + ^27)7 

C27 = 

2(^23 

~ 627 

C28 = 

1(625 + bgg ), 

C29 = 

2(625 

— bgg ). 

C30 = 

2(630 + 631), 

C31 = 

2(^30 

~ 631 


Using the above matrices and relations the computational procedure for calculating Kaluza 
numbers product can be written as follows: 

Y 32 XI = P 32 W 32 A 32 X 512 D 512 P 512 X 32 W 32 P 32 X 32 X 1 (3.3) 

It is easy to see that the diagonal matrix D 512 contains only 32 elements differing by its value. 
The remaining 480 elements coincide with these thirty two elements up to a sign. Let us create a 
column vector C 32 XI = [cq, Ci,..., C 31 ]"’", consisting of the all 32 elements of matrix D 512 which 
possess different values. Then it is easy to see that the elements {ck}, k = 0,1,..., 31 can be 
calculated using the following vectormatrix procedure: 


C32XI — 2 ® H2)B32x1, B32XI — bi, . . . , 631]"'". 
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Figure 2. Data flow diagram for rationalized Kaluza numbers multiplication algorithm in accor¬ 
dance with the procedure (3.3) 
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Fig. 2 shows a data flow diagram representation of the rationalized algorithm for computation 
of Kaluza numbers product and Fig. 3 shows a data flow diagram of the process for calculating the 
vector C 32 XI elements. In Fig. 2, the points where lines converge denote summation. As follows 
from Fig. 3, calculation of elements of diagonal matrix D 512 requires performing only trivial 
multiplications by the power of two. Such operations may be implemented using convention 
arithmetic shift operations, which have simple realization and hence may be neglected during 
computational complexity estimation [l24fl . 


K 

bi 

b?, 



Cl 

C2 

C 3 


^30 C.Q 

^ 31 *^ ^31 


Figure 3. The data flow diagram describing the process of calculating elements of the vector 

C32XI 


4. Evaluation of computational complexity 

We calculate how many real multiplications (excluding multiplications by power of two) and 
real additions are required for realization of the proposed algorithm, and compare it with the 
number of operations required for a direct evaluation of matrixvector product in Eq. (2.3). As 
already mentioned the number of real multiplications required using the proposed algorithm is 
512. Thus using the proposed algorithm the number of real multiplications to calculate Kaluza 
number product is halved. The number of real additions required using our algorithm is 576. 
Therefore, the total number of arithmetic operations for proposed algorithm is approximately 
46% less than that of the direct evaluation. 

5. Conclusions 

In this paper, we have presented an original algorithm allowing to multiply Kaluza num¬ 
bers with reduced both multiplicative and additive complexity. Furthermore, the total number of 
operations in our algorithm is almost two times less than the total number of operations in the 
compared algorithm. Therefore, the proposed algorithm is better than the naive algorithm, both in 
terms of hardware implementation and in terms of its software implementation on a conventional 
computer too. 








21 


References 

[1] Kantor I., Solodovnikov A. Hypercomplex numbers, Springer-Verlag, New York, (1989). 

[2] Billow T., Sommer G. Hypercomplex signals a novel extension of the analytic signal to the multidi¬ 
mensional case, IEEE Trans. Sign. Proc., Vol. SP-49, No 11, 2844-2852, (2001). 

[3] Alfsmann D. On families of 2^-dimensional hypercomplex algebras suitable for digital signal pro¬ 
cessing, in Proc. European Signal Processing Conf. (EUSIPCO 2006), Elorence, Italy, (2006). 

[4] Alfsmann D., Gockler H. G., Sangwine S. J., Ell T. A. Hypercomplex Algebras in Digital Signal 
Processing: Benefits and Drawbacks (Tutorial). Proc. EURASIP 15th European Signal Processing 
Conference (EUSIPCO 2007), Poznan, Poland, 1322-1326, (2007). 

[5] Sangwine S. J., Bihan N. Ee Hypercomplex analytic signals: extension of the analytic signal concept 
to complex signals, Proc. EURASIP 15th European Signal Processing Conference (EUSIPCO 2007), 
Poznan, Poland, 621-624, (2007). 

[6] Moxey C. E., Sangwine S. J., Ell T. A. Hypercomplex correlation techniques for vector images, IEEE 
Trans. Signal Processing, Vol. 51, No 7, 1941-1953, (2003). 

[7] Bayro-Corrochano E. Multi-resolution image analysis using the quaternion wavelet transform. Nu¬ 
merical Algorithms, Vol. 39, No 1-3, 35-55, (2005). 

[8] Calderbank R., Das S., Al-Dhahir N., Diggavi S. Construction and analysis of a new quaternionic 
Space-time code for 4 transmit antennas. Communications in information and systems Vol. 5, No 1, 
1-26, (2005). 

[9] Belfiore J.-C., Rekaya G. Quaternionic lattices for space-time coding. Proceedings of the Information 
Theory Workshop. IEEE, Paris 31 March - 4 April 2003, 267-270, (2003). 

[10] Ertug O. Communication over Hypercomplex Kahler Manifolds: Capacity of Dual-Polarized 
Multidimensional-MIMO Channels. Wireless Personal Communications, vol. 41, No 1, 155-168, 
(2007). 

[11] Silvestrov V. V. Number Systems, Soros Educational Journal, No 8, 121-127, (1998). 

[12] Makarov O. An algorithm for the multiplication of two quaternions, Zh. Vychisl. Mat. Mat. Eiz., Vol. 
17, No 6, 1574-1575, (1977). 

[13] Cariow A., Cariowa G., Algorithm for multiplying two octonions. Radioelectronics and Communi¬ 
cations Systems (Allerton Press, Inc. USA), Vol. 55, Issue 10, 464-473, (2012). 

[14] Cariow A., Cariowa G., An algorithm for fast multiplication of sedenions. Information Processing 
Eetters, 113, 324-331, (2013). 

[15] Cariow A., Cariowa G., ”An algorithm for multiplication of Dirac numbers”. Journal of Theoretical 
and Applied Computer Science, vol. 7, No 4, 26-34, (2013). 

[16] Cariow A., Cariowa G., An Algorithm for Past Multiplication of Pauli Numbers. Advances in Ap¬ 
plied Clifford Algebras, Advances in Applied Clifford Algebras, Vol. 25, No 1, 53-63, (2015). 

[17] Cariow A., Cariowa G., On the Multiplication of Biquaternions, Soft Computing in Computer and 
Information Science: Advances in Intelligent Systems and Computing, Vol. 342, 423^34, (2015). 

[18] Cariow A., Cariowa G., Kubsik B., An algorithm for multiplication of split-octonions, CoRR 
abs/arXiv:1503.01058vl, 1-14, (2015). 

[19] Cariow A., Cariowa G., Knapinski J., Derivation of a low multiplicative complexity algorithm for 
multiplying hyperbolic octonions, CoRR abs/arXiv:1502.06250vl, 1-15, (2015). 



22 


[20] Cariow A., Cariowa G., A new algorithm for multiplying two Dirac numbers, CoRR 
abs/arXiv:1501.00828vl, 1-14, (2015). 

[21] Cariow A., Cariowa G., An algorithm for multiplication of trigintaduonions. Journal of Theoretical 
and Applied Computer Science, Vol. 8, No 1, 50-75, (2014). 

[22] Jariov A., "fariova G., Zracjonalizowany algorytm mnozenia dwoch kwaternionow. Przegl^d Elek- 
trotechniczny. No 9, 137-140, (2010). 

[23] "fariova G., "fariov A., Aspekty algorytmiczne redukcji liczby blokow mnoz^cych w ukladzie do 
obliczania iloczynu dwoch kwaternionow, w czasopismie Pomiary, Automatyka, Kontrola, No 7, 
668-690, (2010). 

[24] "fariov A., "fariova G., ’’Aspekty algorytmiczne organizacji uktadu procesorowego do mnozenia liczb 
Cayleya”, Elektronika, No 11, 137-140, (2010). 

[25] Cariow A. Strategic racjonalizacji obliczeh przy wyznaczaniu iloczynow macierzowo-wektorowych, 
Metody Informatyki Stosowanej, No 1, 147-158, (2008). 

[26] Cariow A. Algorytmiczne aspekty racjonalizacji obliczeh w cyfrowym przetwarzaniu sygnalow, 
Wydawnictwo Zachodniopomorskiego Uniwersytetu Technologicznego, (2011). 



